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We theoretically study the many-body effects of electron electron interaction on the single particle spectral 
function of doped bilayer graphene. Using random phase approximation, we calculate the real and imaginary 
part of the self-energy and hence the spectral function. The spectral function near the Fermi surface shows the 
usual quasiparticle peak, establishing doped bilayer graphene, in contrast to the unstable neutral system, to be 
a Fermi liquid. Away from the Fermi surface, an additional broad plasmaron peak is visible in the spectral 
function. From the low energy behaviour of the self-energy we calculate the quasiparticle residue and the 
effective mass of the quasiparticles as a function of carrier density. We present results for both the on-shell and 
the off-shell approximation for the quasiparticle renormalization. 

PACS numbers: 



Bilayer graphene (BLG) has attracted both experimental 
and theoretical interest as a novel 2D gapless chiral electron- 
hole system with parabolic dispersion (TJ, where carrier den- 
sity can be easily tuned by applying a gate voltage. From a 
many-body interaction physics perspective, BLG is particu- 
larly intriguing as it lies between the chiral and gapless single 
layer graphene (SLG), which has linear band dispersion, and 
2D semiconductor-based electron gas (2DEG) systems, which 
are non-chiral and gapped, but typically have quadratic band 
dispersion [2J. In addition, the interaction parameter in BLG 
can be simply tuned by changing the carrier density whereas 
in SLG the interaction parameter is the fine structure constant 
which is fixed for a given substrate. BLG is therefore the ideal 
system to study chiral interaction physics since interaction ef- 
fects in SLG are typically small [1 1. It is therefore quite sur- 
prising that very little BLG interaction physics has so far been 
studied in the literature, perhaps because of serious techni- 
cal difficulties of treating both intraband and interband virtual 
processes on an equal footing in a gapless system. In this Let- 
ter we provide the first theoretical treatment of the many-body 
interaction effects on the BLG single-particle properties as a 
function of carrier density. Our work, being directly compara- 
ble to experimental data, should motivate experimental studies 
of BLG many-body properties. 

An important many -body property of BLG is the single par- 
ticle spectral function A(k, u), which measures the probabil- 
ity of finding an electron (hole) with a momentum k and en- 
ergy u> in the system (We take h = 1 throughout.) The spectral 
function, which gives detailed information about electronic 
structure and many-body renormalization, can be measured 
by angle resolved photo-emission spectroscopy (ARPES) J3). 
In the absence of interaction, the spectral function is simply a 
delta function giving the non-interacting band energy at mo- 
mentum k. The modification of the spectral function due to 
electron-electron interactions is a matter of great interest both 
theoretically and experimentally. A pressing question in this 
respect is whether BLG behaves like a Fermi liquid, where 
phase space restrictions near the Fermi surface lead to sharply 



defined quasiparticles with a finite spectral weight. 

In this letter, we calculate the carrier density dependent sin- 
gle particle spectral function of doped BLG, taking into ac- 
count electron-electron interactions within the dynamical ran- 
dom phase approximation (RPA). We find that the BLG spec- 
tral function has a broad plasmaron peak, i.e. a peak asso- 
ciated with the quasiparticle carrying a cloud of virtual plas- 
mons, in addition to the Landau quasiparticle peak. We ob- 
tain the quasiparticle renormalization and the effective mass, 
which completely define the low energy coherent part of the 
spectral function around the Fermi surface, and study their 
variation with the carrier density. We find that quasiparti- 
cle spectral weight decreases with decreasing density while 
the effective mass shows a non-monotonic behaviour. The 
effective mass remains smaller than the bare mass down to 
relatively low densities with a maximum renormalization of 
~ 20%, thus showing that the system remains a weakly inter- 
acting Fermi liquid even at low densities. Earlier BLG many- 
body theoretical work restricted itself [4] to either Hartree- 
Fock or static screening theories neglecting dynamical corre- 
lations, which is known to be inadequate in interacting elec- 
tron systems where dynamical screening is important. 

A key question for BLG, indeed for any electronic mate- 
rial, is whether the presence of interaction preserves the Fermi 
liquid behavior, where the system manifests one-to-one corre- 
spondence with the noninteracting system in its low-energy 
properties, or leads to an exotic non Fermi liquid. This ques- 
tion has recently been addressed in the literature for neu- 
tral undoped BLG which exists only at the Dirac point. Not 
surprisingly such a pure zero-density Fermi system is found 
to be a non-Fermi-liquid in the presence of interaction. We 
show in the current work that this Dirac point non-Fermi liq- 
uid behavior of neutral BLG is unstable to the existence of 
any finite doping, and doped BLG is always a Fermi liquid 
in the presence of any finite carrier density (which is always 
true in any physical BLG system). Our finding that the BLG 
is a Fermi liquid similar to regular non-chiral 2DEG systems 
is directly experimentally verifiable through ARPES since the 
singularity structure of the theoretical spectral function is very 
different in the doped and the undoped situations. We find that 
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FIG. 1: (Color online) Left column: The real (thick black line) and 
imaginary part (dashed blue line) of self-energy for (a) k — 0.75/cf, 
(c) k — kF and (e) k = 1.25fc_F- The absolute value of the imaginary 
part is plotted in these graphs. Right column: Single particle spectral 
function for the same wave vectors, (b) k — 0.75/cf, (d) k = kp and 
(f) k — 1.25fcp. The system is characterized by r s — 3. The Dirac 
point is at uj — —Ef- 



the neutral non-Fermi liquid fixed point is an unstable set of 
measure zero, which will always become a Fermi liquid fixed 
point in the presence of density fluctuations. This ia an impor- 
tant new result since density fluctuations are invariably present 
in all real BLG samples due to thermal and disorder effects. 

We work with the low energy two band parabolic approxi- 
mation to the BLG Hamiltonian where the dispersions of the 
bands are given by E sk — sk 2 /2rriQ where s = ±1 de- 
notes the conduction and valence bands and mo is the non- 
interacting BLG mass with mo — 0.033m e , m e being the 
mass of a free electron. The chiral band wave functions are 
given by ip\ v = (e l2e *, s)/V2, where 6 = tan" 1 (k y / k x ) 
is the azimuthal angle in the momentum space. Each band 
has a degeneracy factor of g = 4 (2 for spin and 2 for val- 
ley). We note that the 2-band parabolic BLG approximation 
is well-justified because we are interested in the low-energy 
behavior, and at high doping densities, where the parabolic 
approximation fails, the BLG dispersion becomes linear with 
the system behaving like the SLG whose many-body proper- 
ties have already been studied in the literature |[T]|6). 

The non-interacting Green function for the system is given 
by G (s,k,ui) = (u> - &k)"\ where £ sk = E sk - 
Ep with Ep being the Fermi energy. The correspond- 
ing spectral function consists of delta function peaks at the 



band dispersion energies and is given by Aq(s,]s., uj) = 
— (l/7r)JmGo(s, k, w) = 5(u> — £ s k)- Electron-electron in- 
teractions cause scattering of the non-interacting excitations 
resulting in the broadening of the delta function peaks and 
shifting of spectral weight. This is incorporated through the 
modification of the interacting Green function by a self en- 
ergy term, G _1 (s,k,w) = G " 1 (s,k,w) — E s (k,w) and the 
corresponding spectral function is given by 



s;'(k >W ) 



u-s;(k, w )p + [s;'(k, w )p 



(i) 



where S s (k,w + i0+) = E a (k,u;) + i£' s '(k,w). Thus the 
imaginary part of self-energy incorporates the broadening of 
the spectral function peak while the real part determines the 
shift in the energy dispersion, especially near the Fermi sur- 
face. 

We calculate the self-energy and hence the spectral function 
of a doped BLG system within the dynamical RPA at zero 
temperature. Using Matsubara frequencies, the self-energy 
can be written as 



G (s', k - q, iui - iu')F ss , (k, q) 



(2) 



where ioj and iu>' are respectively fermionic and bosonic Mat- 
subara frequencies. Here V c (q) — 2ire 2 /{nq) is the 2D bare 
Coulomb potential (k being the background dielectric con- 
stant). Using Ep as units of energy and fcj? as units of mo- 
mentum, the strength of the interaction can be written in terms 
of the dimensionless coupling parameter r s = e 2 'gin/ '(nkp). 
Unlike SLG (but similar to 2DEG), r s in BLG can be tuned 
by changing the density of carriers (r s ~ ?i -1 / 2 ) and we will 
present our results as a function of r s rather than the density. 
The dielectric function e(q, uj), which incorporates the effects 
of dynamic screening, was calculated analytically in Ref. 7 
Finally, the chirality of the bands is encoded in the overlap 
vertex factor F ss r (k, q) = (1 + ss' cos 26) /2, where 8 is the 
angle between k and k + q. 

Although it is customary to separate the self-energy into 
an exchange and a correlation contribution, we prefer to work 
with the full expression since the exchange contribution, by it- 
self, is divergent due to presence of the filled Fermi sea in the 
valence band. The total contribution, however, is free of this 
divergence, showing that the use of the dynamically screened 
interaction is crucial in calculating the single particle spec- 
tral properties of BLG. It is important to emphasize that our 
dynamical RPA self-energy calculation is theoretically well- 
defined and does not involve any arbitrary infrared or ultravi- 
olet regularization through arbitrary momentum cut-offs. 

Spectral function: With no loss of generality, we assume 
the BLG chemical potential to be in the conduction band. 
Since we are interested in the low energy properties around 
the Fermi surface, we will focus on A(+, k, uj) and drop the 
index s from now on. In the left hand column of Fig. [T]we 
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FIG. 2: Single particle spectral function for the wave vectors, (a) 
k — 0.75fcF, (b) k — kF and (c) k = 1.25&F for a system character- 
ized by r s = 7. The low energy quasiparticle peaks are accompanied 
by plasmaron features away from k = fcp. 



plot the calculated real and imaginary parts of the self-energy 
for three wavevectors, (a) k = 0.75kp (b)fc = fop and (c) 
k = 1.25kp, in a system characterized by r s — 3, corre- 
sponding to carrier density n = 3x 10 12 cm -2 for S1O2 sub- 
strate. The self-energy has two distinct observable structures, 
coming from scattering off plasmon collective modes (plas- 
marons) and the continuum of particle hole excitations. In the 
imaginary part of the self-energy, the plasmon contribution 
leads to distinct logarithmic singularities away from k = kp. 
These singularities are accompanied by jumps in the real part 
of the self-energy. The continuum contribution leads to an in- 
creasing background and dominates the self-energy at large 
values of a;. In the right hand column of Fig. [TJ we plot the 
corresponding BLG spectral function. We also show the spec- 
tral function for a larger value of r s — 7 (lower density of 
n = 6 x 10 n cTO~ 2 ) in Fig. [2] The spectral functions at these 
two values of r s show qualitatively similar features, although 
for the larger value of r s , the incoherent background part of 
the spectral weight is spread over a larger energy. The promi- 
nent feature of the spectral function, manifesting its Fermi 
liquid characteristics, consists of a quasiparticle peak which 
disperses across the Fermi energy as the wave-vector goes 
across the Fermi surface. The quasiparticle peak is broad- 
ened for wavevectors away from k = kp. At k = kp, the 
spectral function has a delta function peak at the Fermi en- 
ergy with a weight Z ~ 0.4 for r s — 3 and Z ~ 0.28 for 
r s = 7. This is the classic Fermi liquid behavior of an inter- 
acting electron system. Away from the Fermi wavevector, the 
spectral function consists of a double peaked structure. The 
second peak, arising from the dressing of the quasiparticles 
by the plasmon modes, is often referred to in the literature 
as a plasmaron. In SLG systems, these satellite bands have 
already been observed in ARPES experiments |3), and our 
current work shows that both the quasiparticle peak and the 



plasmaron should be manifestly observable in BLG as well. 

Landau's great insight in to the theory of interacting Fermi 
systems was that due to phase space constraints arising from 
Pauli blocking, the imaginary part of the self energy van- 
ishes rapidly as one approaches the Fermi surface and weakly 
interacting quasiparticles form the low energy elementary 
excitations of the system. In a 2DEG, at low energies, 
the width of the quasiparticle peak Tk = ImJ^(k, £k) ~ 
(£k/-^'i r ) l n (£k/-E'F)- This ensures that the spectral function 
has a delta function peak at the Fermi surface with a finite 
weight Z < 1. For BLG, we analytically find that the leading 
order quasiparticle scattering rate follows the same asymptotic 
form, as shown in Fig. [3ja). In fact this can be analytically 
shown to be true for arbitrarily large values of r s . Thus the 
system remains a Fermi liquid for arbitrarily small densities. 
Several authors [5 | have recently shown that the strictly un- 
doped zero density neutral graphene is a non Fermi liquid with 
a vanishing Z factor and 7mS(fc, uj) ~ uj. However this non 
Fermi liquid behaviour relies on two peculiar features of un- 
doped BLG : (a) a quadratic dispersion around a single Fermi 
point and (b) the interband scattering contribution persisting 
upto uj — 0. Both these conditions are violated at any fi- 
nite carrier density: the dispersion becomes linear around the 
Fermi surface and the interband scattering contribution is cut- 
off at a scale of Ep. The low energy quasiparticle width is in 
fact dominated by forward scattering from intraband particle- 
hole continuum and follows the usual Fermi liquid form. In 
the renormalization group language, the chemical potential is 
a relevant perturbation (even at the tree level) and thus the non 
Fermi liquid fixed point is unstable to arbitrarily small densi- 
ties. We therefore conclude that the BLG is always a Fermi 
liquid. 

Thus, for the doped BLG near the Fermi surface, the spec- 
tral function can be written as 

A(k, uj) = Z k S(uj - v F k + + A mc (k, uj) (3) 

where /1 is the chemical potential of the system, vp is the 
slope of the dispersion, Z is the quasiparticle spectral weight 
and A lnc the non-singular incoherent part of the spectral func- 
tion. The effective mass of the quasiparticles, which enters 
into low temperature thermodynamic properties like specific 
heat, is defined as m* — kp/vp. The self-energy at low ener- 
gies thus contains important information about the quasiparti- 
cle dispersion, broadening and spectral weight. 

Quasiparticle spectral weight: As mentioned before, phase 
space restrictions associated with the fermionic Pauli principle 
lead to sharply defined quasiparticles near the Fermi surface in 
the interacting system producing the one-to-one Fermi liquid 
correspondence with the non-interacting Fermi gas. In BLG, 
the imaginary part of the self-energy at k = kp vanishes as 
uj 2 In uj as one approaches the Fermi energy and exactly at the 
Fermi surface, the spectral function has a delta function peak 
with a weight Z < 1. This leads to a discontinuity in the 
momentum distribution at the Fermi surface (at T = 0) with 
the magnitude of the discontinuity given by Z. In terms of the 
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FIG. 3: (Color online): (a) Width of the quasiparticle peaks, Tk = 
JmE(k, £k), as a function of k/kF for two different values of r a . 
The red dashed ine is for r s — 3 and the thick black line is for 
r s = 7. Inset: The low energy part of Tk/^k as a function of £k 
for r s = 3 and r s — 7. The x-axis is logarithmic and the linear 
plot shows the lj 2 In uj scaling of the quasiparticle width, (b) The 
quasiparticle weight, Z, as a function of r s . (b) The effective mass 
m* as a function of r a . In (b) the thick black line is obtained from 
the on-shell approximation whereas the blue dashed line is the self- 
consistent solution. 



self energy the quasiparticle spectral weight is given by 



z- x = \- 



dT, (fcjr, a>), 



\k F ,E^ 



(4) 



In Fig. [3|a), we plot the BLG quasiparticle weight as a func- 
tion of r s . The quasiparticle weight decreases monotoni- 
cally with increasing r s (decreasing density), as stronger in- 
teractions shift spectral weight away from the quasiparticle 
pole. For very large r s , Z eventually slowly approaches zero 
asymptotically with the system remaining a Fermi liquid for 
all r s except the zero density infinite r s singular point. 

The energy dispersion of the quasiparticle excitations 
can be extracted either from the self-consistent solution of 
Dyson's equation (often called the off-shell approximation) 



Wfc - £k - £ (k,w fe ) = 
or from the on-shell approximation (OSA) 

c^-£ k -£'(k,&) = 



(5) 



(6) 



It is clear that if the self-energy is calculated exactly, the 
Dyson equation must be solved self-consistently. However 
the RPA involves a single-loop leading-order approximation 



in the dynamically screened Coulomb interaction and as such, 
the dispersion should be calculated within the OSA to get con- 
sistent results within RPA [8]. However, since the use of the 
off-shell approximation is found in various places in the liter- 
ature, we present our results both within the OSA and the off- 
shell Dyson equation approximations. We find that the inter- 
action effects for the off-shell approximation are smaller than 
those for the OSA. We believe that experimental measure- 
ments of effective mass and/or the Z-factor compared with 
our theory could shed light on the important question of the 
relative validity of on-shell versus off-shell approximation. 

The effective mass (to*): The slope of the quasiparticle dis- 
persion around the Fermi surface determines the Fermi veloc- 
ity vf which can be written in terms of the effective mass of 
the Landau quasiparticles as vp = kp/m* . The self consis- 
tent solution to the Dyson equation yields 



TO 

to 



Z~ 



1 



OS (k.i 



(7) 
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where Z is the quasiparticle weight. The on-shell approxima- 
tion gives 
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1 + k F dk \k F ,E F i- 



(8) 
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The effective mass obtained from both calculations is plotted 
as a function of r s in Fig. |3jb). We find that the effective 
mass at first decreases with increasing r s and then exhibits a 
weak nonmonotonicity for larger r s , where it tends to increase 
a bit, with the renormalized mass being always less than the 
bare band mass even for r s as large as 10. This is very differ- 
ent from 2DEG and 3DEG systems where the effective mass 
renormalizes to values much larger than the bare mass for r s 
larger than unity |8|. This prediction of our theory, which is 
somewhat counterintuitive (i.e. the renormalized mass being 
less than the bare mass even for strong interaction), is eas- 
ily experimentally verifiable. We note that the BLG effective 
mass renormalization is small (only of the order of ~ 20% 
or less) in sharp contrast to 2DEG or 3DEG where the mass 
renormalization could be by a factor of 5 or more at larger 
r s J8). The relatively small BLG mass renormalization indi- 
cates a larger regime of validity for dynamical RPA in the chi- 
ral gapless BLG system compared with the standard 2DEG 
system. The maximum mass renormalization of ~ 20% also 
shows that the system remains a weakly interacting Fermi liq- 
uid upto a relatively large r s . 

Conclusion: We have calculated the effects of electron- 
electron interaction on the single particle spectral function 
of bilayer graphene within the random phase approximation. 
The spectral function shows a clear satellite plasmaron peak 
in addition to the Landau quasiparticle peak near the Fermi 
surface. Thus, it should be possible to obtain information 
about collective plasmon excitations in the system by look- 
ing for these satellite peaks in ARPES data. Focusing on the 
low energy properties of the Landau quasiparticles, we have 



5 



calculated the quasiparticle weight and the effective mass of 
the Landau quasiparticles as a function of r s . We show that 
inclusion of dynamic screening is crucial in obtaining these re- 
sults. We also show that the presence of the valence band leads 
to relatively small renormalization in the low energy spectral 
properties, and thus the system remains a weakly interacting 
Fermi liquid upto relatively large values of r s or relatively 
low densities. In addition to obtaining the quasiparticle spec- 
tral function, the effective mass renormalization, the inelastic 
scattering rate, and the quasiparticle renormalization factor, 
we have established the theoretical principle that doped BLG 
is always a Fermi liquid in contrast to undoped BLG. 
This work is supported by ONR-MURI and NRI-SWAN. 
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